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Abstract We evaluate the measurement precision of the pro¬ 
duction cross section times the branching ratio of the Higgs 
boson decaying into tau lepton pairs at the International Lin¬ 
ear Collider (ILC). We analyze various final states associ¬ 
ated with the main production mechanisms of the Higgs 
boson, the Higgs-strahlung and WW-fusion processes. The 
statistical precision of the production cross section times 
the branching ratio is estimated to be 2.6% and 6.9% for 
the Higgs-strahlung and WW-fusion processes, respectively, 
with the nominal integrated luminosities assumed in the ILC 
Technical Design Report; the precision improves to 1.0% 
and 3.4% with the running scenario including possible lu¬ 
minosity upgrades. The study provides a reference perfor¬ 
mance of the ILC for future phenomenological analyses. 

1 Introduction 

After the discovery of the Higgs boson by the ATLAS and 
the CMS experiments at the LHC [1,2], the investigation 
of the properties of the Higgs boson has become an im¬ 
portant target of study in particle physics. In the Standard 
Model (SM), the coupling of the Higgs boson to the matter 
fermions, i.e., the Yukawa couplings, is proportional to the 
fermion mass. The Yukawa couplings can deviate from the 
SM prediction in the presence of new physics beyond the 
SM. Recent studies indicate that the deviations from the SM 
could be at the few-percent level if there is new physics at 
the scale of around 1 TeV [3]. It is therefore desired to mea¬ 
sure the Higgs couplings as precisely as possible in order to 
probe new physics. 

In this study, we focus on Higgs boson decays into tau 
lepton pairs (h —t+t^) at the International Linear Collider 
(ILC). This decay has been studied by the ATLAS and the 
CMS experiments, who reported a combined signal yield 
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consistent with the SM expectation, with a combined ob¬ 
served significance at the level of 5.5(7 [4-6]. The purpose 
of this study is to estimate the projected ILC capabilities of 
measuring the decay mode in final states resulting 

from the main Higgs boson production mechanisms in 
collisions. Existing studies on h -A decays at 
collisions [7,8] did not take into account some of the rele¬ 
vant background processes or were based on a Higgs boson 
mass hypothesis which differs from the observed value, both 
of which are addressed in this study. We assume the ILC ca¬ 
pabilities for the accelerator and the detector as documented 
in the ILC Technical Design Report (TDR) [9-13] together 
with its running scenario published recently [14, 15]. The 
results presented in this paper will be useful for future phe¬ 
nomenological studies. 

The contents of this paper are organized as follows. In 
Section 2, we describe the ILC and the ILD detector con¬ 
cept, and the analysis setup. The event reconstruction and 
selection at center-of-mass energies of 250 GeV and 500 
GeV are discussed in Sections 3 and 4. Section 5 describes 
the prospects for improving the measurement precision with 
various ILC running scenarios. We summarize our results in 
Section 6. 


2 Analysis conditions 

2.1 International Linear Collider 

The ILC is a next-generation electron-positron linear col¬ 
lider. It covers a center-of-mass energy {^/s) in the range of 
250-500 GeV and can be extended to — 1 TeV. In the 
ILC design, both the electron and positron beams can be po¬ 
larized, which allow precise measurements of the properties 
of the electroweak interaction. The details of the machine 
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Table 1 Typical integrated luminosities L and center-of-mass energies 
v/iofthelLC [9,14,15], 


Scenario 

v/7 (GeV) 

L(fb-') 

Nominal 

250 

250 


500 

500 

Luminosity upgrade 

250 

2000 


500 

4000 




Fig. 1 Diagrams of the main production mechanisms of the Higgs bo¬ 
son in collisions. Top: Higgs-strahlung {ZK) process. Middle: 

VFW-fusion process. Bottom: ZZ-fusion process. 


design are summarized in the ILC Technical Design Report 
(TDR) [9-13], 

The ILC aims to explore physics beyond the SM via pre¬ 
cise measurements of the Higgs boson and the top quark as 
well as to search for new particles within its energy reach. 
The center-of-mass energies and integrated luminosities which 
are foreseen are summarized in Table[T] The numbers for the 
nominal running scenario are taken from the ILC TDR [9]. 
The numbers for scenarios including energy and luminosity 
upgrades are based on studies in Refs. [14,15]. 


2.2 Production and decay of the Higgs boson 

Figure [T] shows the diagrams for the main production mech¬ 
anisms of the Higgs boson in collisions. The cross sec¬ 
tions of Higgs boson production calculated by WHIZARD [16] 
with a Higgs mass of 125 GeV are shown in Figure|2] where 
polarizations of—80% and -1-30% for the electron and positron 
beams are assumed, and initial state radiation is taken into 
account. 

For the calculation of the production cross section and 
the subsequent decay of the signal processes of ^ 



\/s (GeV) 


Fig. 2 Cross sections of the Higgs boson production as a function of 
s/s in the electron-positron interaction. 


ffh —^ where / denotes a fermion, we use an event 

generator based on GRACE [17,18]. The effect of beamstrahlung 
is implemented according to the calculation by GuineaPig [19], 
which simulates beam-beam interactions, with the beam 
parameters described in the TDR [20]. Initial state radiation 
is incorporated following the prescription developed by the 
ILC Event Generator Working Group [ 13,21 ]. To handle the 
spin correlation of tau pairs from the Higgs boson decay, 
GRACE is interfaced with TAUOLA [22-25]. The decays of 
other short-lived particles and the hadronization of quarks 
and gluons are handled by PYTHIA [26]. 


2.3 Background processes 

For background processes, we use common Monte-Carlo 
(MC) samples for SM processes previously prepared for the 
studies presented in the ILC TDR [13]. The event samples 
include > 2/, -> 4/, > 6 /, and ^ 

ffh. The event generation of these processes is performed 
with WHIZARD [16], in which beamstrahlung, initial state ra¬ 
diation, decay of short-lived particles, and hadronization are 
taken into account in the same way as described in the previ¬ 
ous section for the signal process. The background processes 
from 77 interactions with hadronic final states, in which pho¬ 
tons are produced by beam-beam interactions, are gener¬ 
ated on the basis of the cross section model in Ref. [27]. We 
find that the interactions between electron or positron beams 
and beamstrahlung photons, i.e., e* 7 —e* 7 , e* 7 —?> 3/, and 
5/, have negligible contributions to background. 
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Table 2 Summary of the performance of the ILD detector model. 


Name 

Value 

Impact parameter resolution 

c 10 

5© - 57 T— pm 

psin^©0 


, 1 X 10-^ 

Momentum resolution 

Jet energy resolution 

2 X 10 5 ^ GeV/c 

Pt Sin 6 

30 

-7^ 

^E{Ge<) 


2.4 Detector Model 

The detector model used in this analysis is the International 
Large Detector (ILD), which is one of the two detector con¬ 
cepts described in the ILC TDR. It is a general-purpose An 
detector designed for particle flow analysifl aiming at best 
possible jet energy resolution. 

The ILD model consists of layers of sub-detectors sur¬ 
rounding the interaction point. One finds, from the inner¬ 
most to the outer layers, a vertex detector (VTX), a silicon 
inner tracker (SIT), a time projection chamber (TPC), a sili¬ 
con envelope tracker (SET), an electromagnetic calorimeter 
(ECAL), and a hadron calorimeter (HCAL), all of which are 
put inside a solenoidal magnet providing a magnetic field 
of 3.5 T. The return yoke of the solenoidal magnet has a 
built-in muon system. The ILD design has not yet been fi¬ 
nalized. In this analysis, we assume the following configu¬ 
rations and performance. The VTX consists of three double 
layers of silicon pixel detectors with radii at 1.6 cm, 3.7 cm 
and 6 cm. Each silicon pixel layer provides a point resolu¬ 
tion of 2.8 /fm. The TPC provides up to 224 points per track 
over a tracking volume with inner and outer radii of 0.33 m 
and 1.8 m. The SIT and SET are used to improve the track 
momentum resolution by adding precise position measure¬ 
ments just inside and outside of TPC. The ECAL consists 
of layers of tungsten absorbers interleaved with silicon lay¬ 
ers segmented into 5x5 mm^ cells, has an inner radius of 
1.8 m, and has a total thickness of 20 cm corresponding to 
24 radiation length. The HCAL consists of layers of steel 
absorbers interleaved with scintillator layers segmented into 
3x3 cm^ cells and has an outer radius of 3.4 m correspond¬ 
ing to 6 interaction length. Additional silicon trackers and 
calorimeters are located in the forward region to assure her¬ 
metic coverage down to 5 mrad from the beam line. The 
key detector performance of the ILD model is summarized 
in Table |2] Details of the ILD model and the particle flow 
algorithm are found in Refs [13,28]. 

*The particle flow algorithm aims at achieving the best attainable jet 
energy resolution by making one-to-one matching of charged particle 
tracks with calorimetric clusters so as to restrict the use of calorimetric 
information, which is in general less precise than tracker information, 
to neutral particles. This requires highly granular calorimeters and a 
tracking system with high performance pattern recognition for events 
with high particle multiplicity. 


2.5 Detector simulation and event reconstruction 

In this study, we assume a Higgs boson mass of 125 GeV, 
a branching ratio of the Higgs boson decay into tau pairs 
(BR(/! —T+T^)) of 6.32% [29], and beam polarizations of 
—80% and 4-30% for the electron and the positron beams, 
respectively. 

We perform a detector simulation with Mokka [30], a 
Geant4-based [31] full detector simulator, with the ILD model 
for all signal and background processes, with the excep¬ 
tion of the 4- 2/ process at -Js = 500 GeV, 

for which SGV fast simulation [32] is used. The event re¬ 
construction and physics analysis are performed within the 
MARLIN software framework [33], in which events are re¬ 
constructed using track finding and fitting algorithms, fol¬ 
lowed by a particle flow analysis using the PandoraPFA 
package [28]. 

3 Analysis at the center-of-mass energy of 250 GeV 

At -y/i = 250 GeV, the Higgs-strahlung (e+e^ — Zh) pro¬ 
cess dominates the SM Higgs production, as shown in Eig- 
ure|2l The VTVT-fusion and ZZ-fusion cross sections are neg¬ 
ligible at this energy. We take into account ffh (ex¬ 
cluding the A —!► signal), ^ 2 /, and ^ 

4/ for the background estimation. The 77 —^ hadrons back¬ 
ground is overlaid onto the MC samples with an average of 
0.4 events per bunch crossing [27]. An integrated luminosity 
of 250 fb^^ is assumed for the results in this section. 

There are four main signal modes: qqh, 

e+e^A, /i+/i^A, and vvA. Eor our ^ = 

250 GeV results, we report on the first three of these modes. 
We do not quote the results for the vvA mode as we And that 
it suffers from background processes with neutrinos in the 
final state. We do not analyze the t+t^A mode in 

this study. 

3.1 —>■ qqh 

Reconstruction of isolated tau leptons and the Z ^ qq 
decay 

Eor the qqh mode, we first identify the tau leptons us¬ 
ing a dedicated algorithm developed for this topology. The 
algorithm proceeds as follows. 

1. The charged particle with the highest energy is chosen 
as a working tau candidate. 

2. The tau candidate is combined with the most energetic 
particle (charged or neutral) satisfying the following two 
conditions: the angle 0 , between the particle and the tau 
candidate satisfies cos 0 , > 0.99; and the combined mass, 
calculated from the sum of the four momenta of the par¬ 
ticle and the tau candidate, does not exceed 2 GeV. The 
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Fig. 3 Distributions of the invariant mass of the reconstructed tau lep¬ 
ton pairs at y/s = 250 GeV for the -5- qqh mode. M^+^- and Mcoi 
stand for the tau pair masses before and after the collinear approxima¬ 
tion, respectively, for the signal. M(;oi(bkg) is the tau pair mass with the 
collinear approximation for the background. 


four momentum of this particle is then added to that of 
the tau candidate. 

3. Step 2 is repeated until there are no more particles left to 
combine. The resulting tau candidate is then set aside. 

4. The algorithm is repeated from Step 1 until there are no 
more charged particles left. 

A tau candidate is accepted if the number of charged parti¬ 
cles with track energy greater than 2 GeV is equal to one or 
three, the net charge is equal to ± 1 , and the total energy 
is greater than 3 GeV. Furthermore, an isolation require¬ 
ment is applied as follows. A cone of half-angle 0^, with 
cos 0e = 0.95, is dehned around the direction of the tau mo¬ 
mentum. The tau candidate is accepted if the energy sum of 
all particles inside the cone (excluding those forming the tau 
candidate) does not exceed 10 % of the tau candidate energy. 
We require exactly two hnal tau candidates with opposite 
charges. This results in a selection efficiency of 49.3% for 
the signal events. 

After the tau candidates are identihed, the neutrino en¬ 
ergy is recovered by using the collinear approximation [34]. 
Because tau leptons from a Higgs boson decay are highly 
boosted, it is reasonable to assume that the tau momentum 
and the neutrino momentum are nearly parallel. Under this 
assumption, the energy of the two neutrinos, one from each 
tau decay, can be solved by requiring that the overall trans¬ 
verse momentum of the event is balanced in two orthogonal 
directions. The neutrino reconstructed in this way is added to 
the tau candidate. Figure [3 shows the invariant mass distri¬ 
butions of the tau pairs without and with (Mcoi) the 

collinear approximation for the events containing two tau 
lepton candidates with opposite charges. With the collinear 
approximation, a clear peak is visible at 125 GeV for signal 
events. The M^oi distribution for background events with the 
same criteria is also shown. 


The Durham jet clustering algorithm [35] is applied to 
the remaining particles to reconstruct the two jets from the 
Z boson decay. 

Event selection 

We perform a pre-selection over the reconstructed events, 
followed by a multivariate analysis. The pre-selection is de¬ 
signed to reduce background while keeping most of the sig¬ 
nal. The events are pre-selected according to the following 
criteria. The Z ^ qq candidate and the h —^ candidate 

are successfully reconstructed. The total number of charged 
particles is at least 9. The visible energy of the event, Eyis, 
lies in the range of 105 GeV < Eyis < 255 GeV. The visi¬ 
ble mass of the event, Myis, is greater than 95 GeV. The sum 
of the magnitude of the transverse momentum of all visible 
particles, Pf,sum, is greater than 40 GeV. The thrust of the 
event is less than 0.97. The Z candidate dijet has an energy, 
Ez, in the range of 60 GeV < Ez < 175 GeV and has an in¬ 
variant mass, Mz, in the range of 35 GeV < Mz < 160 GeV. 
The angle between the two jets, Ojj, satishes cos Ojj < 0.5. 
The recoil mass against the Z boson, computed as Mrecoii = 
y' {^/s - Ezf - |pzp, is in the range of 65 GeV < Mrecoii < 
185 GeV. The Higgs candidate tau pair before the collinear 
approximation has an energy, Et+t ’ 140 GeV and 

an invariant mass, M^+^- , in the range of 5 GeV < < 

125 GeV. The angle between the two tau candidates, 0^+^-, 
satishes cos 9^+^- < —0.1. The tau pair after the collinear 
approximation has an energy, Ecob in the range of 30 GeV 
< Ecoi < 270 GeV and an invariant mass, Mcoi, in the range 
of 15 GeV < Mcoi < 240 GeV. 

We use a multivariate analysis using Boosted Decision 
Trees (BDTs) as implemented in the Toolkit for Multivariate 
Data Analysis [36] of the ROOT framework [37]. The input 
variables are 

- Eyis, Pf.vis, cos 0miss, where P; yis is the magnitude of the 
visible transverse momentum and Omiss is the angle of 
the missing momentum with respect to the beam axis; 

- Mz, cos 9jj, Mi-ecoib cos 0z, where 9z is the angle of the 
Z candidate momentum with respect to the beam axis; 

- M^+t--, cos0^+t.-, cos0acop, where 0acop is the 

acoplanarity angle between the two tau candidates; 

- E logioMo/cUol, E logioko/oiol> where t/o/o-rfo ™tl 

T+,T“ T+,T“ 

zo/cUo respectively the transverse and longitudinal 
impact parameters of the most energetic track in the tau 
candidate divided by their respective uncertainty esti¬ 
mated from the track ht; 

- Mcoi, Ecob cos0cob where 0coi is the angle of the 
Higgs candidate momentum with the collinear approxi¬ 
mation measured from the beam axis. 

The BDTs are trained using a set of statistically independent 
signal and background samples. The distribution of the re¬ 
sulting multivariate discriminant is shown in Figure |4] We 
apply a hnal selection on the multivariate discriminant that 
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Fig. 4 Distributions of the multivariate discriminant from training 
Boosted Decision Trees for the > qqh mode, shown for the sig¬ 

nal and the total background. 

Table 3 Event yields estimated for the -5- qqh mode at y/s = 250 
GeV, assuming an integrated luminosity of 250 fb“* and beam po¬ 
larizations of P(e“,e+) = (—0.8,-1-0.3), shown for the signal and the 
background processes. The signal contribution {h —> T+T^) is removed 
from the ffh process. “No cut” is the number of events corresponding 
to the production cross section times the integrated luminosity. “Pre¬ 
selected” is the number of events after the pre-selection for the multi¬ 
variate analysis. “Final” is the number of events after the selection on 
the multivariate discriminant. 



Signal 

ffh 

2/ 

4/ 

No cut 

3318 

7.649 X 10^ 

2.863 X 10’ 

1.736 X 10* 

Pre-selected 

1451 

3526 

2316 

6.940 X lO'^ 

Final 

1232 

22.0 

9.3 

512.0 


maximizes the signal significance defined as 5 /^yS + B, where 
S and B are the number of signal and background events, 
respectively. The hnal selected sample consists of 1232 sig¬ 
nal and 543 background events. The estimated event yields 
before and after the selection are summarized in Table [3] 
The signal selection efficiency is 37% with a signal signifi¬ 
cance of 29, which corresponds to a statistical precision of 
A{ax BR )/{a x BR) = 3.4%. 

3.2 —>■ e^e^h 

Z boson and tan lepton reconstruction 

For the e^e^h mode, we hrst reconstruct the pair 
that forms a Z boson candidate. A reconstructed particle is 
identified as an electron or a positron if its track momen¬ 
tum (Ptrk) and its associated energy deposits in the ECAL 
(Eecal) and HCAL (Ehcal) satisfy the following criteria; 

Eecal/(Eecal TEhcal) > 0.96, 

(Thecal + £’hcal)/A rk > 0.6. 

For the particles that are identified as electrons or positrons, 
we further require that \dQ /| < 6 and |zo/cEo I < 3, to re¬ 
duce the electrons from secondary decays such as the tau 


lepton decays from the Higgs boson. We also require the 
track energy to be greater than 10 GeV, which removes the 
contamination from the 77 —> hadron background. The e+e— 
pair whose combined mass is closest to the Z boson mass is 
selected as the Z boson candidate. To improve the mass and 
energy resolutions, the momenta of nearby neutral particles 
are added to that of the Z candidate if their angle 0 mea¬ 
sured from at least one of the satisfies cos 0 > 0.999. 
The fraction of signal events that survive the Z 

boson selection is 61%. 

We apply a tau finding algorithm to the remaining par¬ 
ticles. Compared with the qqh mode, the algorithm is sim¬ 
pler due to the absence of hadronic jet activities aside from 
the tau decays. Starting with the charged particle with the 
highest energy as a working tau candidate, we define a cone 
around its momentum vector with a half-angle of 0^ = 1.0 
rad. Particles inside the cone are combined with the tau can¬ 
didate if the combined mass remains smaller than 2 GeV. 
The tau candidate is then set aside, and the tau hnding is re¬ 
peated until there are no more charged particles left. The tau 
candidates are then separated into two categories according 
to its charge. Within each category, the tau candidate with 
the highest energy is selected. The chosen pair forms 
the Higgs candidate. Finally, the collinear approximation is 
applied to the selected tau candidates. 

Event selection 

A pre-selection is applied with the following require¬ 
ments before proceeding with the multivariate analysis. The 
Z —candidate and the h candidate are suc¬ 

cessfully reconstructed. The total number of charged tracks 
is 8 or fewer, which ensures statistical independence from 
the q'qh mode. The visible energy is in the range of 100 
GeV < Fvis < 280 GeV. The visible mass is in the range 
of 85 GeV < Mvis < 275 GeV. The sum of the magnitude 
of the transverse momentum of all visible particles, A.sum, is 
greater than 35 GeV. The Z —candidate has an energy 
in the range of 40 GeV < Ez < 160 GeV and an invariant 
mass in the range of 10 GeV < Mz < 145 GeV. The recoil 
mass against the Z boson, Mfecoii, is greater than 50 GeV. 

We then apply a multivariate analysis using BDTs using 
the following input variables: 

- Mvis, Fvis, cos 0miss, COS 0thrust, where 0thmst is the angle 

of the thrust axis with respect to the beam axis; 

- Mz, Mfecoiri 

- cos 0T+r-, COS 0acop; 

- L logioMo/OdJ, and £ logio ko/c^zol- 

T+,T“ 

A final selection on the multivariate discriminant is applied 
to maximize the signal significance, giving 76.3 signal and 
44 background events. The final signal selection efficiency 
is 44%. The estimated event yields before and after the se¬ 
lection are summarized in Tabled The signal significance is 
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Table 4 Event yields estimated for the mode at y(s = 

250 GeV, assuming an integrated luminosity of 250 fb“* and beam 
polarizations of P(e^,e+) = (—0.8,+0.3). Refer to Tablej^for the row 
definitions. 



Signal 

ffh 

2/ 

4/ 

No cut 

175.1 

7.964 X 10“* 

2.863 X 10’ 

1.736 X 10* 

Pre-selected 

109.4 

60.2 

3.334 X 10^ 

1.169 X lO'^ 

Final 

76.3 

4.2 

0 

39.9 


estimated to be 7.0, corresponding to a statistical precision 
of 4 ((7 X BR)/((7 X BR) = 14.4%. 

3.3 ^ 

Z boson and tan lepton reconstruction 

The reconstruction procedure of this mode is similar to 
that of the e^e^h mode, with the electron identification re¬ 
placed by the muon identification. The muons are identified 
by requiring 

Thecal/(^ecal + £’hcal) < 0.5, 

(^ECAL +£’HCAL)/Rtrk < 0.6. 

We additionally require the identified muons to satisfy \do /o^g | 
3 and |zo/<7zol ^ ^ track energy greater than 

20 GeV. The efficiency for selecting such muon pairs in 
/f^ T+ signal events is 92%. The tau lepton reconstruc¬ 
tion is the same as in the e^e^h mode. 

Event selection 

The following pre-selection requirements are applied be¬ 
fore proceeding with the multivariate analysis. TheZ — 
candidate and the are successfully reconstructed. 

The total number of charged tracks is 8 or fewer. The vis¬ 
ible energy is in the range of 105 GeV < Evis < 280 GeV. 
The visible mass is in the range of 85 GeV < Myis < 275 
GeV. The sum of the magnitude of the transverse momen¬ 
tum of all visible particles, is greater than 35 GeV. 

The Z ^ candidate has an energy in the range of 45 

GeV < Ez < 145 GeV and an invariant mass in the range 
of 25 GeV < Mz <125 GeV. The recoil mass against the 
Z boson, Mfecoib is greater than 75 GeV. The invariant mass 
of the tau pair system before the collinear approximation, 
is smaller than 170 GeV. 

A multivariate analysis with BDTs is applied to the pre¬ 
selected events using the following input variables: 

“ 4^vis5 Eyis, 7^?,yis5 

- Mz, COS dz, iWrecoii; 

- cos Mqq\, 

- L logiol'^o/oi/gl, and ^ logio |zo/oVo|. 

We apply a final selection on the multivariate discriminant 
to maximize the signal significance, and obtain 101.9 signal 


Table 5 Event yields estimated for the -5- mode at y/s = 

250 GeV, assuming an integrated luminosity of 250 fb“* and beam 
polarizations of P(e^,e+) = (—0.8, +0.3). Refer to Table[3]for the row 
definitions. 



Signal 

ffh 

2/ 

4/ 

No cut 

164.6 

7.965 X 10“* 

2.863 X 10’ 

1.736 X 10* 

Pre-selected 

132.8 

63.5 

4182 

8011 

Final 

101.9 

2.2 

0 

29.0 


and 31 background events. The final signal selection effi¬ 
ciency is 62%. The estimated event yields before and after 
the event selection are shown in Table |5] The signal signif¬ 
icance is estimated to be 8 . 8 , corresponding to a statistical 
precision of A (<7 x BR) /(o’ x BR) = 11.3%. 

4 Analysis at the center-of-mass energy of 500 GeV 

At i/s = 500 GeV, both the WW-fusion and the Higgs-strahlung 
processes have sizable contributions to the total signal cross 
section. We take into account the —>■ ffh (except h 

T+T^), —>■ 2/, —>■ 4/, and 6 / processes 

as backgrounds. The 77 -+> hadron background is overlaid 
onto the signal and background MC samples, assuming an 
<average rate of 1.7 events per bunch crossing [27]. The anal¬ 
ysis in this section assumes an integrated luminosity of 500 
fb^'. We report our results on the qqh and 

vvh modes. We do not give results for the e^e^h 

and modes, as they do not contribute sig¬ 

nificantly to the overall sensitivity due to their small cross 
sections. 

4.1 —>■ qqh 

Reconstruction of isolated tau leptons and the Z ^ qq 
decay 

We start with the tau finding, following the same proce¬ 
dure described in Section 3.1. We additionally require the 
tau candidate to have an energy greater than 4 GeV. The en¬ 
ergy of the neutrino from tau decays is corrected using the 
collinear approximation as before, resulting in a clear peak 
around the Higgs boson mass as can be seen in Figure|5] We 
find that 54% of qqT^T^ signal events survive the require¬ 
ment of finding exactly one pair of T+ . 

The invariant mass of all the particles, except those be¬ 
longing to the two identified tau candidates, should be con¬ 
sistent with the Z boson mass; however, a shift to a higher- 
mass value is observed, due to the presence of non-negligible 
background particles from 77 —hadron events contaminat¬ 
ing signal events. In order to mitigate the effect of these 
background particles, we use the kr clustering algorithm [38, 
39] implemented in the Fast Jet package [40] with a gen¬ 
eralized jet radius of R = 0.9. The jets that are formed along 
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e*e'^qqh,\ls=500GeV,L=500fb\P(e',e*)=(-80%,+30%) 



Fig. 5 Distributions of the invariant mass of the reconstructed tau lep¬ 
ton pairs at y/s = 500 GeV for the -5- qqh mode. M^+^- and Mcoi 
stand for the tau pair masses before and after the collinear approxima¬ 
tion, respectively, for the signal process. Mcoi(bkg) is the tau pair mass 
with the collinear approximation for the background processes. 

the beam axis are then discarded. The remaining particles 
are clustered into two jets by using the Durham clustering 
algorithm to reconstruct the Z boson decay. 

Event selection 

To facilitate the multivariate analysis, we impose the fol¬ 
lowing pre-selections. The Z ^ qq candidate and the h —?> 
candidate are successfully reconstructed. The total num¬ 
ber of charged tracks is between 8 and 70. The visible energy 
of the event is in the range of 140 GeV < Eyis < 580 GeV. 
The visible mass of the event is in the range of 120 GeV 

< Mvis < 575 GeV. The sum of the magnitude of the trans¬ 
verse momentum of all visible particles, Pf,sum, is greater 
than 70 GeV. The thrust of the event is less than 0.98. The Z 
candidate dijet has an energy in the range of 50 GeV <Ez < 
380 GeV and has an invariant mass in the range of 5 GeV 

< Mz < 350 GeV. The recoil mass against the Z boson is 

in the range of 40 GeV < Mrecoii < 430 GeV. The Higgs 
candidate tau pair before the collinear approximation has an 
energy, l^ss than 270 GeV and an invariant mass, 

less than 180 GeV, and the angle between the two 
tau candidates satisfies cos0^+^- < 0.7. The tau pair after 
the collinear approximation has an energy in the range of 40 
GeV < Ecoi < 430 GeV and an invariant mass, Mcoi, which 
is less than 280 GeV. 

A multivariate analysis with BDTs is applied using the 
following input variables: 

- Evis^ ^r,sum, ^vis^ where Pvis is the magnitude of the visi¬ 
ble momentum; 

- Mz, Ez, cos Ojj, cos Mrecoii; 

- M^+^-, COS 0^+^-, Mcoi, 

- E logic Mo/OrfJ, and E logic ko/oij. 

After choosing the optimum threshold on the multivariate 
discriminant to maximize the signal significance, we are left 
with 782 signal and 335 background events. The final sig¬ 


nal selection efficiency is 37%. The event yields before and 
after the selection are summarized in Table|6l The signal sig¬ 
nificance is found to be 23.4, corresponding to a statistical 
precision of A (<7 x BR) /(o’ x BR) = 4.3%. 

4.2 —>■ vvh 

Tau pair reconstruction 

The tau finding algorithm proceeds in the same way as 
described for the e^e^h mode in Section 3.2, except that the 
half-angle of the cone Qc around the most energetic track 
is modified to 0.76 rad. The most energetic positively and 
negatively charged tau candidates are combined to form a 
Higgs boson candidate. 

Event selection 

For the vvh mode, it is necessary to suppress the large 
background coming from the -f 2/ processes. 

We apply the following requirements to mitigate this back¬ 
ground. A tau lepton pair T+ is successfully reconstructed. 
The total number of tracks is less than 10. There is at least 
one charged track with a transverse momentum greater than 
3 GeV and at least one charged track with an energy greater 
than 5 GeV. The missing momentum angle with respect to 
the beam axis satisfies |cos0miss| < 0.98. The acoplanarity 
angle between the two tau candidates satisfies cos 0acop < 
0.98. At this point, 94% of the +2/ back¬ 

ground is eliminated, while retaining 85% of the signal events 

The following additional pre-selections are applied be¬ 
fore the multivariate analysis. The visible energy is in the 
range of 10 GeV < Eyis < 265 GeV. The visible mass is in 
the range of 5 GeV < Myis < 235 GeV. The missing mass, 
4^missj is greater than 135 GeV. The sum of the magnitude 
of the transverse momentum of all visible particles, 
is greater than 10 GeV. The Higgs candidate tau pair be¬ 
fore the collinear approximation has an energy, less 

than 240 GeV and an invariant mass, of less than 

130 GeV. The angle between the two tau candidates satisfies 
cos 0T+T- < 0.8. A requirement on the transverse impact pa¬ 
rameter of the tau candidate which gives a smaller value of 
the two is applied, such that min|^fo/( 77 (,| > 0.01. 

A multivariate analysis with BDTs is applied using the 
following input variables: 

- Number of tracks with energy greater than 5 GeV; 

- Number of tracks with transverse momentum greater than 

5 GeV; 

“ 4^ViS5 Eyis, yis, Ef sumi cos 0thrust, COS ^miss5 

- COS COS 0acop5 

- logic min Ic/o/o-rfJ. 

We obtain 1642 signal and 1.11 x 10^ background events af¬ 
ter optimizing the selection on the multivariate discriminant. 
The final signal selection efficiency is 30%. The event yields 
before and after the selection are summarized in Table |7] 














Table 6 Event yields estimated for the e+e —> qqh mode at ^Js = 500 GeV, assuming an integrated luminosity of 500 fb * and beam polarizations 
of P{e^ ,e^) = (—0.8, +0.3). Refer to Table[3]for the row definitions. 


_ Signal ffh _2/_4/_ 6/ 

No cut 2131 1.266 xlO^ 1.320x 10’ 9.989x 10* 6.929 x 10= 

Pre-selected 1088 2889 3.013x 10“’ l.UdxlO^ 1.737x 10’ 

Final 782.1 17.6 1.5 275 41 


Table 7 Event yields estimated for the e+e —5- v v/t mode at y/s = 500 GeV, assuming an integrated luminosity of 500 fb ’ and beam polarizations 
of P{e^ ,e^) = (—0.8, +0.3). Refer to Table[3]for the row definitions. 


_ Signal ffh _2/_4/_ 6/ 

No cut 5534 1.232x 10^ 1.320x 10’ 9.989x 10* 6.929 x 10^ 

Pre-selected 3623 1543 5.957x 10“* 1.756x 10’ 990.8 

Final 1642 65.5 379 1.043x 10’ 238 


e*e-^vVh,\/s=500GeV,L=500fb',P(6',e*)r(-80%,+30%) 



Fig. 6 Distribution of neutrino-pair invariant mass, Myy, using event 
generator information for the > vv/i mode at y/s = 500 GeV. 

The effect of detector simulation is not applied. 


The signal significance is 14.5, corresponding to a statistical 
precision of 4 (a x BR)/((7 x BR) = 6.9%. 

In this mode, the Zh ^ vvh process and the 

—7> vvh process via ITlT-fusion are expected to be 
the dominant contributions. The effect of the interference 
between these two processes is studied using the distribu¬ 
tion of the invariant mass of the neutrino pair Myv com¬ 
puted from event generator information as shown in Fig¬ 
ure |6l A clear peak around the Z boson mass is visible, 
with a small contribution underneath it coming from the tail 
from higher masses, indicating that the interference of the 
Zh —^ vvh process and the ITVT-fusion is small. 
We hence split the events into two categories based on this 
generator-level variable, and define events with Myv < 120 
GeV as Z/z” events, and those with Myv >120 

GeV as “VTVT-fusion” events. We find that the 1642 signal 
events after the final selection is composed of 13% —?> 

Zh events and 87% VTW-fusion events. The selection effi¬ 
ciencies for Zh and WVT-fusion events are 33.5% 

and 29.2%, respectively. 


5 Discussion 

5.1 Precision with the ILC running scenarios 

We now discuss the prospects of the measurement precision 
with the ILC running scenarios proposed in Refs. [14, 15] 
by extrapolating the results presented in the previous sec¬ 
tions. Table |8] summarizes the integrated luminosities for 
various center-of-mass energies and beam polarizations for 
three different scenarios we consider. 

In order to estimate the statistical precision of the cross 
section times the branching ratio measurements with elec¬ 
tron and positron beam polarizations other than (—0.8, +0.3) 
used in the previous sections, we need to know the cor¬ 
responding selection efficiencies for the signal and back¬ 
ground processes. In the following, we assume the same 
selection efficiencies obtained in the previous sections for 
all of these beam polarizations, although in principle the 
angular distributions of the final states may depend on the 
beam polarizations. This assumption is nevertheless justi¬ 
fied as follows. The Zh process is mediated by the 

s-channel Z boson exchange with the vector or the axial 
vector coupling, which forbids the same-sign helicity states 
(± 1, ± 1), while giving more or less the same angular distri¬ 
butions for the opposite-sign helicity states (=F 1, ± 1). On the 
other hand, the VTVT-fusion process proceeds only through 
the left-right helicity states (—1,+1), since the W boson 
couples only to the left-handed and the right-handed e+. 

For the signal processes, therefore, their angular distribu¬ 
tions stay the same for the active (i.e. opposite) helicity states, 
independently of the choice of beam polarizations. The same 
reasoning applies to the background processes with the s- 
channel yjZ exchange or those involving W bosons coupled 
to the initial state e+ or . On the other hand, the processes 
involving f-channel photon exchange or photon-photon in¬ 
teractions do not forbid the same-sign helicity states. How¬ 
ever, since the probability of finding an electron and a positron 
in the same-sign helicity states is the same for both the (—0.8, +0.3) 
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Table 8 Integrated luminosity (in fb *) for various beam polarizations (P[e ), P(e+)) at = 250 GeV and 500 GeV for three running scenarios. 
“Nominal” is the scenario described in the TDR [9] and used in the main results of this paper. “Initial” and “Full” are the scenarios proposed in 
Refs. [14,15]. 


Scenario 

vA (GeV) 

(-80,+30) 

(+80, -30) 

(-80,-30) 

(+80,+30) 

Total 

Nominal 

250 

250 

0 

0 

0 

250 


500 

500 

0 

0 

0 

500 

Initial 

250 

337.5 

112.5 

25 

25 

500 


500 

200 

200 

50 

50 

500 

Full 

250 

1350 

450 

100 

100 

2000 


500 

1600 

1600 

400 

400 

4000 


and (+0.8, —0.3) beam polarizations, the efficiencies for such 
background processes with the same-sign helicity states should 
also be the same. In our estimation, we do not use the re¬ 
sults of the (±0.8, ±0.3) beam polarizations, since the sig¬ 
nal cross sections are small and the integrated luminosities 
collected at these beam polarizations are foreseen to be small. 
Under these assumptions, the selection efficiencies will not 
depend on the choice of beam polarizations. We can then es¬ 
timate the projected statistical precision for other scenarios 
by calculating the number of signal and background events 
with the production cross sections and the integrated lumi¬ 
nosities for individual beam polarizations, according to the 
running scenarios. The result from this estimation is sum¬ 
marized in Tabled 

5.2 Precision of the h —branching ratio 

So far, we discussed the precision of the production cross 
section times the branching ratio, which is the primary in¬ 
formation we will obtain from the experiments. Here, we 
discuss the prospects for measuring the branching ratio it¬ 
self. At the ILC, the production cross section for the Higgs- 
strahlung process can be separately measured using the re¬ 
coil mass technique [10,13]. The cross section for the WW- 
fusion process can also be determined by using the branch¬ 
ing ratio for the h^bb decay [10]. The obtained cross sec¬ 
tion values allow us to derive the branching ratio for the 
—?► decay. 

At ^/s = 250 GeV, the contributions of the WW-fusion 
and ZZ-fusion processes are negligible. Therefore, we can 
use the Higgs-strahlung cross section to derive the branch¬ 
ing ratio. The Higgs-strahlung cross section Ozh can be mea¬ 
sured to a statistical precision of Aozhl^^zh = 2.5% with the 
nominal TDR running scenario [13]. This improves to a sub¬ 
percent level with the full running scenario [14]. 

At i/s = 500 GeV, both the Higgs-strahlung and the VTVT- 
fusion processes contribute to the Higgs boson production, 
whereas the contribution of the ZZ-fusion process is negli¬ 
gible. For the ^ vvh mode, in which both processes 

are present, it is in principle possible to estimate the contri¬ 
butions from the Higgs-strahlung and the VTVT-fusion pro¬ 


cesses separately, as discussed in Section 4.2. However, we 
do not use this mode here for the estimate. The expected sta¬ 
tistical precision of the branching ratio after combining all 
the modes except the vvfi mode is 3.6% for the nominal run¬ 
ning scenario. This improves to 1.4% with the full running 
scenario, where we assume Aozh/’^zh = 1-0%. 

5.3 Systematic uncertainties 

The MC statistical uncertainties are found to have negligi¬ 
ble impact on the results. The systematic uncertainty in the 
luminosity measurement has been estimated to be 0.1% or 
better for the ILC [41] and is not expected to be a significant 
source of systematic errors. The uncertainties in the selec¬ 
tion criteria, such as those caused by the uncertainty in the 
momentum/energy resolutions and tracking efficiencies, are 
not included in this analysis, since they are beyond the scope 
of this paper. 

6 Summary 

We have evaluated the measurement precision of the Higgs 
boson production cross section times the branching ratio of 
decay into tau leptons at the ILC. The study is based on the 
full detector simulation of the ILD model. The dominant 
Higgs boson production mechanisms were studied at the 
center-of-mass energies of 250 GeV and 500 GeV, assum¬ 
ing the nominal luminosity scenario presented in the ILC 
TDR. The analysis results are then scaled up to the running 
scenarios taking into account realistic running periods and a 
possible luminosity upgrade. 

The results for the various modes and scenarios are sum¬ 
marized in Table |9] In short, the cross section times the 
branching ratio can be measured with a statistical precision 
of A (cT X BR)/(<7 X BR) = 2.6% and 1.0% for the nominal 
and full running scenarios, respectively. We evaluate the sta¬ 
tistical precision of BR(/i —^ T+T^) to be 3.6% for the nom¬ 
inal TDR integrated luminosity and 1.4% for the full run¬ 
ning scenario, respectively. These results serve to provide 
primary information on the expected precision of measuring 
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Table 9 Expected precision of the cross section times the branching ratio 4((7 x BR)/(c7 x BR), assuming various running scenarios. 


Scenario 

v/J (GeV) 

L(fb-') 

qqh 

e^e h 

h 

vvh 

Combined 

Nominal 

A(ffxBR)/((TxBR) 

250 

250 

3.4% 

14.4% 

11.3% 

— 

3.2% 


500 

Combined 

Combined 

500 

4.3% 

2.7% 

14.4% 

11.3% 

6.9% 

6.9% 

2.6% 

6.9% 

Initial 

A(ffxBR)/((TxBR) 

250 

500 

2.5% 

10.9% 

8.7% 

— 

2.4% 


500 

Combined 

Combined 

500 

4.9% 

2.3% 

10.9% 

8.7% 

9.6% 

9.6% 

2.1% 

9.6% 

Pull 

A(ffxBR)/((TxBR) 

250 

2000 

1.3% 

5.5% 

4.3% 

— 

1.2% 


500 

combine 

combine 

4000 

1.7% 

1.0% 

5.5% 

4.3% 

3.4% 

3.4% 

1.0% 

3.4% 


Higgs decays to tau leptons at the ILC, which will be useful 
for future phenomenological studies on physics beyond the 


SM. 
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